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We study numerically the finite temperature and frequency mobility of a 
particle coupled by a local interaction to a system of spinless fermions in one 
dimension. We find that when the model is integrable (particle mass equal 
to the mass of fermions) the static mobility diverges. Further, an enhanced 
mobility is observed over a finite parameter range away from the integrable 
point. We present a novel analysis of the finite temperature static mobility 
based on a random matrix theory description of the many-body Hamiltonian. 
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The description of dissipation in many-body systems is a problem that traditionally 
attracted the interest of theorists [[[]]. In particular, in classical systems, a relation was 
observed between the occurence of normal heat transport and chaotic behavior [0] . Recently 
we have found a similar relation in the finite temperature transport properties of a quantum 
many-body system [[5]. This system, describing a particle coupled to a fermion environment 
@], shows dissipationless transport at finite temperatures when integrable while normal 
transport when non-integrable. 

The analysis of the particle mobility was based on the Kubo linear response theory ||; 
within this framework, an indication of dissipationless transport is given by a non-vanishing 
charge stiffness Q D = lim w ^ wfi"(w)/2, where /i" is the imaginary part of the mobility. We 
have established that, at any finite temperature T > 0, D remains non-zero when the model 
is integrable, while D decreases exponentially to zero with system size in the generic non- 
integrable case 0. We attributed the non- vanishing D to integrability motivated by recent 
studies relating transport properties to level statistics and level statistics to integrability 
[T0| , [TT| . In particular, D can be expressed as the thermal average of level curvatures of the 



system subject to an external flux and thus it reflects the properties of the energy spectrum. 

However integrability, and the resulting non-vanishing charge stiffness, is a singular effect 
in parameter space. In this work, we are going to study if dissipationless transport is signaled 
by an enhanced static mobility in the neighborhood of the integrable point. In order to 
promote an effect as a signature of an integrable point, its robustness to small perturbations 
is indeed an important issue. 

The model we are studying, describes the evolution of a single particle on a chain of L 
sites with periodic boundary conditions. The particle is coupled to an environment of iV 
spinless fermions through a local repulsive interaction. The Hamiltonian is : 

H = -t h J2(d] +1 di + H.c.) -tJ2(c] +1 ci + H.c.) + (1) 
i i i 

where q(cJ) are annihilation (creation) operators for the fermions at site I and di(d\) for the 



particle. The model is integrable by Bethe's ansatz for th = t and non-integrable for 



th 7^ t; we have verified that the level-spacing distribution is Poisson for the integrable case 
and follows the Gaussian Orthogonal Ensemble (GOE) distribution for the non-integrable 
one. 

From linear response theory, the real part of the mobility fi'(w) is : 

1 — 

n'(u) = 2nD8(u) +7T Pn \ (n\j \ m)\ 2 5(uj - e m + e n ). (2) 

where | n), e n denote the eigenstates and eigenvalues of the Hamiltonian, p n = exp(—/3e n )/Z 
the Boltzmann factors, (3 = 1/T, Z the partition function and j = ith J2i (d] +1 di — d]di + \j 
the current that does not commute with H. The f-sum rule relates the mobility to the average 
kinetic energy K of the particle : /f^ fj,(u))duj = tt(—K) where K = —thYliW+idl + H.c.) 

Before dealing with the fine effect of integrability, we can comment on the general fre- 
quency behavior of the mobility. For large interaction parameter U, fi'(w) consists of two 
well separated peaks at energies u ~ and to ~ U. At high temperatures, the weight and 
the width of the high energy peak can be estimated from moments to be ir(—K) p(l — p) and 
yz(3t| + Qt 2 ) respectively |T^]; the weight of the peak is essentially given by the number 
of particle-hole pairs in the fermion bath while its width depends on the hopping matrix 
elements. For smaller U, the two peaks merge into a single peak at low frequencies. 

Below we will present a study of fi'(uj) from a numerical simulation for finite systems. For 
the smallest sytems, (L < 12), the energies and matrix elements of the current are calculated 
exactly using a full-diagonalization technique. For larger systems, (L = 14,16), fi'(uj) is 
evaluated using a method proposed recently for finite temperature dynamical quantities, 
combining random sampling of states with Lanczos diagonalization |15j . In most cases, we 
will present results for the integrated mobility I(u) = fi'(u)')du)' /(2tc(— K)), lo > 0, thus 
avoiding the usual artificial broadening of the 5-functions in Eq.(0). I (to) has two simple 
limits, lim^o — D/ (-K) and lim^oo I(uo) = 1/2. The slope of gives the mobility 
itself. 

Our investigation is restricted to the high temperature regime, to relatively small U = 2t 
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and at half filling, p = N/L = 1/2. The high temperature regime is indeed easier to 
analyse numerically since the finite size effects are less pronounced and the choice of U allows 
to maximize the effect of integrability since the charge stiffness decreases with increasing 
U. In the following we present numerical results for n'(u>), discuss the static mobility, 
fi Q = lim^o ti>'{w) around the integrable point and estimate the finite size effects. We then 
propose a new approach for evaluating /io based on a recent relation between dissipation and 
random matrix theory (RMT) f|||. We conclude with a presentation of /io as a function of 
t h - 

In Fig. we present /j,'(uj) for different size systems at T — 5t and far from the integrable 
point. We note that fi'(u>) is well defined with no excessive size dependence (the restricted 
number of Lanczos iterations in the approximate method does not allow us to reach the 
lowest frequencies for L = 14,16). The size dependence is a very important issue since 
normal transport can be established only if /io is intensive, i.e. if it saturates with system 
size 0. The range of available sizes is rather limited but nevertheless we observe that /i 
has mostly saturated for L — 12. Therefore we can expect to have real transport in the 
high temperature regime even for such small systems. The indicated full points at u = 0, in 
this and the following figures, are the estimation of /i by the RMT approach that will be 
discussed below. 

Fig. H presents I{oS) for L = 12 and T — 5t scanning th through the integrable point 
th = t. We first note that at high frequencies, u > 1.5t, the particularity of the integrable 
point is not reflected in I(oS). On the other hand, we observe that, at low frequencies, all 
the weight is concentrated in the charge stiffness for th = t, while it spreads over a finite 
frequency range for th^t. 

We now focus on the static mobility, around th = t. It is a very difficult quantity to 
extract from numerical results on finite systems because of the discreteness of the spectrum. 
We found convenient to estimate /io from the uj —>■ limit of (I(ou) — 1(0)) /uj which is 
presented in Fig. |3| for different th- We note that the low frequency limit is well defined for 
th > 1.5t] it becomes more difficult to extract /io for th — t as (I(u>) — l(0))/u varies rapidly 
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close to ijj = 0. Nevertheless, we see that the extrapolated — 1(0)) /u gives a value of 

Ho that clearly increases as t^ approaches the integrable point t^ = t. 

In contrast, we would like to illustrate the low frequency behaviour of for th — t. 

Fig. ^ shows I{oS) for different system sizes up to L = 16. The stars at u — indicate the 
exact charge stiffness in the thermodynamic limit as was computed analytically in a previous 
work using the Bethe-ansatz method ||. The weight at low frequencies is concentrated in 
the charge stiffness I(uo = 0) and reaches this point with a zero slope indicative of a 
fjLo = 0. Although we still observe a residual size dependence, we think that the depletion 
of weight at low frequencies persists in the thermodynamic limit since I(u — 0) is already 
close to the analytical value. When U is decreased the shoulder at uj ~ t moves to lower 
frequencies and the depletion of weight is observed in a smaller frequency range. These 
observations suggest that the effect of integrability shows up in the low frequency range, i.e. 
for uj lower than all the natural energy scales, t, U. 

From a first point of view we would expect the effect of integrability, absence of level 
repulsion and modification of level density - density correlations, to be reflected only in a 
frequency range of the order of the level spacing A (in which case it might be amenable 
to a RMT analysis as in e.g. [jl6|| ). However, in our many-body system, A scales to an 
exponentially small value with the size of the system (of the order of 10~ 2 t for L = 12). The 
reason that our results for //(o>) show the influence of integrability at u ^> A, is the strong 
dependence of the current matrix elements on \th — 1\ in this frequency range. Nevertheless, 
as we will show in the next section, we can still obtain an estimate of /xq by assuming that, 
locally in the spectrum, the current matrix is described by a random matrix. 

Ho from random matrix theory : In the non-integrable case, we expect the Hamiltonian 
matrix of our system to be described, at least locally, by random matrices from the GOE 
ensemble (as we verified for the level- spacing distribution). We can then rewrite the Kubo 
formula for the static mobility /io using a relation between off-diagonal and diagonal elements 
of the current for random matrices in this ensemble || : 
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2 (I Jnm | 2 We m = ((inn " v{E)f) = 0*{Ef) (3) 

Here j nm — (n \ j | m), v = (j nn ) and the outer brackets denote an average of matrix elements 
locally in the spectrum, i.e. for e n ,e m ~ E (within a given translational symmetry sector). 
Numerical tests show that indeed, for th ^ t, the above relation (|3|) is well satisfied. From 
(0) the static mobility p$ can be expressed in terms of the variance a 2 and the many-body 
density of states p(E) at the energy E : 

(a sum over all translational symmetry sectors is implied). This expression may be obtained 
from (I(5(jj) —I(0))/8u, for 5uj ^> A but much smaller than the natural energy scales, which 
is the appropriate limit for finding the static response in a macroscopic system. 

Further, as the diagonal elements of the current can be calculated through the first 
derivatives of the energies with respect to a fictitious flux <f>, j nn = de n /d<j), we see that (|J) 
relates a transport property, po, to a change in boundary conditions; <ft is introduced in the 
Hamiltonian by the usual Peierls construction th — > the 1 ^. Similar relations have been used 



in the theory of localization JL7[ and in the study of the conductivity of metallic systems 
at zero temperature (see also [18]). We propose it here for a many-body system at finite 
temperature . 

We use expression (|j) to evaluate the static mobility for our model in the high tem- 
perature regime, (3 — > 0. Our finite systems do not allow a test of this expression at low 
temperatures because there might appear deviations from GOE statistics in the low energy 
part of the spectrum. The po so extracted was shown in Fig. [I] and Fig. [3] by the full 
points at uo = in fair agreement with the extrapolated value from p'{uj). These results give 
support to the RMT approach to transport for our model. In Fig. [5] we collect our results 
for po as a function of th for different size systems. For a fixed size, we observe a divergence 
of po at th — t which becomes more pronounced with increasing size. This indicates that 
in the infinite size limit an enhanced po appears over a finite parameter range around the 
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integrable point. Further study of the critical behaviour at th = t would require larger size 
systems. 

In consistence with a finite po, the above analysis can also give some insight into the 
suppression of the charge stiffness for our system as follows : the second derivative of the 
free energy, F, with respect to the flux is : 

d 2 F 



In this expression the third term vanishes by symmetry and we find that in a system as 
ours, without persistent currents, vanishes in the thermodynamic limit. Therefore the 
charge stiffness is also given by : 
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If the current matrix can be described by a random matrix locally in the spectrum, the 
variance a 2 has to scale with system size as l/p(E); this is consistent with a finite po f° r 
the thermodynamic limit in Eq.(^j). Moreover we observed that v is microscopic in the bulk 
of the spectrum. Therefore, the charge stiffness itself is exponentially suppressed in the 
thermodynamic limit system since p(E) oc exp(aL). It would be interesting to see whether 
the RMT approach for the finite temperature static conductivity also applies to other generic 
non-integrable systems in the metallic regime. 

In conclusion, the numerical simulations presented in this work allowed us to establish 
a scenario for the signature of integrability in the finite frequency mobility. In general, the 
dependence of p'iu) is governed by U, th and t; at low frequencies, the weight is concentrated 
in a peak centered at zero frequency of roughly constant magnitude but of decreasing width 
as the integrable point is approached. For = t this peak collapses to a 5-function with 
weight equal to the charge stiffness. As the parameter range over which an enhanced static 
mobility is observed is finite, this study suggests that the effect of integrability on transport 
is a physically relevant effect. 
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FIGURES 

FIG. 1. Normalized finite frequency mobility for th = 2t, U = 2t, T = 5t and different system 
sizes. The full point at to = indicate the static mobility from the RMT approach for L = 14 and 
/9-»0. 

FIG. 2. Integrated mobility I(u) for different t h sX U = 2t, T = 5t and L = 12 : t h = t (thick 

solid line), 0.5t (solid line), 0.7t (dotted line), 1.2t ( dashed line ) and 2i (long dashed line). 

FIG. 3. (/(w) — 1(0)) /uj as a function of frequencies for L = 12, U = 2t, T = 5t and several 

non-integrable cases in the low frequency regime. The full points at uj = indicate the static 

mobility from the RMT approach. 

FIG. 4. Integrated mobility I(u>) for the integrable case th = t at U = 2t and different temper- 
atures and sizes : L = 16 ( solid line ), L = 14 ( dotted line ), L = 12 ( dashed line ) and L = 10 ( 

long dashed line ). The stars indicate the exact charge stiffness in the thermodynamic limit. 

FIG. 5. Normalized static mobility as a function of t^/t in the high temperature limit and for 

U = 2t. The values are computed within the RMT approach. 
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